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ABSTRACT 

We present the results of a numerical investigation of the spectra and light curves of the emission from the 
photospheres of long-duration gamma-ray burst jets. We confirm that the photospheric emission has high ef- 
ficiency and we show that the efficiency increases slightly with the off-axis angle. We show that the peak 
frequency of the observed spectrum is proportional to the square root of the photosphere's luminosity, in agree- 
ment with the Amati relation. However, a quantitative comparison reveals that the thermal peak frequency 
is too small for the corresponding total luminosity. As a consequence, the radiation must be out of thermal 
equilibrium with the baryons in order to reproduce the observations. Finally, we show that the spectrum inte- 
grated over the emitting surface is virtually indistinguishable from a Planck law, and therefore an additional 
mechanism has to be identified to explain the non-thermal behavior of the observed spectra at both high and 
low frequencies. 

Subject headings: Gamma-ray burst: general — radiation mechanisms: thermal — methods: numerical — 
relativistic processes 



1. INTRODUCTION 

Despite more than thirty years of investigation, the origin 
of the prompt emission of GRBs is still veiled in mystery. 
Among the various proposed radiation mechanisms are syn- 
chrotron emission either in a baryon- or magnetically domi- 
nated jet (Piran 1999; Lloyd & Petrosian 2000; Zhang & Yan 
2010), synchrotron self-Compton emission (Pe'er & Waxman 
2004; Baring & Braby 2004), quasi-thermal comptonization 
(Ghisellini & Celotti 1999), bulk Compton emission (Lazzati 
et al. 2000), and photospheric emission, i.e., the radiation 
advected by the outflow that is released as the flow becomes 
optically thin (Goodman 1986; Rees & Meszaros 2005; Peer 
et al. 2005, 2006; Giannios 2006; Lazzati et al. 2009). 

In recent years, the photospheric model has gained consen- 
sus for its robust modeling, its lacking of any adjustable pa- 
rameters, its high efficiency (Lazzati et al. 2009, hereafter pa- 
per I; Mizuta et al. 2010; Nagakura et al. 2010) and its ability 
to easily reproduce the observed peak frequency (Pe'er et al. 
2005. 2006; paper I). Theoretical investigation has moreover 
revealed that in the case of sub-photospheric dissipation, pho- 
tospheric radiation is characterized by prominent non-thermal 
high-frequency tails, due to inverse Compton (IC) scattering 
of the thermal photons off relativistic electrons (Pe'er et al. 
2006; Giannios 2006; Giannios & Spruit 2007; Lazzati & 
Begelman 2010). However, due to the relatively small size 
of their emission region, photospheric models cannot explain 
the very high frequency emission observed by Fermi (Abdo 
et al. 2009; Zhang et al. 2010). GeV photons cannot be 
produced in photospheric models and their origin has to be 
identified with the early afterglow (or external shock) emis- 
sion (Kumar & Barniol Duran 2009; GhiselUni et al. 2010). 
In addition, photospheric models cannot naturally reproduce 
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the non-thermal low-frequency power-law spectra below the 
peak (see, e.g., Pe'er & Ryde 2010). 

In this paper we present an extension of the simulations 
and results of Paper I. We have performed a simulation in a 
10 times larger box, extending to 2.5 x lO'-' cm, allowing us 
to explore the properties of the off-axis emission and to per- 
form a quantitative comparison of our light curves and spec- 
tra with observational data. We also address the issue of the 
low-frequency spectrum which is observed to be non-thermal 
(Kaneko et al. 2006) but is predicted, in the simple thermal- 
photosphere model, to be thermal. In particular we explore 
whether the superposition of spectra of different temperatures 
over the emitting surface can significantly alter the thermal 
appearance of the spectrum (Pe'er & Ryde 2010). 

This paper is organized as follows: in Sect. 2 we describe 
our simulations, in Sect. 3 we show the light curves and spec- 
tra extracted from our simulations, in Sect. 4 we compare our 
synthetic light curves and spectra to the Amati relation, and 
in Sect. 5 we discuss our results in the context of theoretical 
predictions and simulations from other groups. 

2. NUMERICAL SIMULATIONS 

The numerical simulation on which this paper is based is 
an extension of the one presented in paper I. A jet with lu- 
minosity L = 5.33 X 10^" erg/s is injected in the core of a 16 
solar mass Wolf-Rayet star evolved to pre-explosion (model 
16TI, Woosley & Heger 2006). The central engine has a con- 
stant luminosity for 100 s, after which it is sharply turned off. 
The jet is injected at a distance of 10^ cm from the star cen- 
ter with an opening angle Oq = 10°, a Lorentz factor Fq = 5, 
and enough internal energy to reach an asymptotic Lorentz 
factor Poo = 400, upon complete non-dissipative acceleration. 
The jet evolution was computed with the special-relativistic, 
adaptive mesh refinement code FLASH (Fryxell et al. 2000), 
as modified by the authors (Morsony et al. 2007) in a box 
of size 2.6 x lO'"* cm in the jet direction and 5 x 10'^ cm in 
the equatorial direction. Science frames were extracted every 
0.2 seconds. The physics included, the resolution, and the re- 
finement schemes were analogous to those of Morsony et al. 
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Fig. 1 . — Photospheric light curves for different viewing angles. Thick lines show Ught curves obtained directly from the simulations, while thin lines show 
light curves when a correction for the finite size of the simulations box is applied (see text). 



(2007) and paper I. 



3. PHOTOSPHERIC LIGHT CURVES 

The location of the photosphere in the observer direction 
(which we define as the z axis) Zpi, was calculated analo- 
gously to paper I by back-integrating the optical depth in 
space and time from a virtual observer located at a distance 
of Zobs = 1.3 X lO" cm. For each time step in the observer 
frame, the path of a photon is traced back in time as the photon 
travels in the evolving density and velocity patterns. As a re- 
sult the location of the photosphere is computed as a function 
of the observed time, and no transformation from laboratory 
to observed coordinates is required. Note also that we use a 
Cartesian coordinate z for the photosphere rather than the ra- 
dial distance r in order to automatically take into account the 
curvature effect and the equal arriving time surfaces. Let x 
and y be Cartesian coordinates in the plane perpendicular to 
the line of sight and z be the Cartesian coordinate along the 
line of sight. We define x to be the coordinate in the simu- 
lation plane and y the coordinate perpendicular to the simu- 
lation plane (and to the line of sight). Given the calculated 
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comoving density distribution n'(ti^\,,x,z), the location of the 
photosphere for the photons observed at the time fobs and at a 
location x is given by: 



1 



CTT-n fobs 



,x,z] T[l-l3cos(9,)]dz 



(1) 

where (3 = /3{t\ab,x,z) is the local velocity of the outflow in 
units of the speed of light, F = T(t\.^\,,x,z) is the local bulk 
Lorentz factor, and 0,, = 9t.{ti^b,x,z) is the angle between the 
velocity vector and the direction of the line of sight. All the 
values of /3, F, and 6y are evaluated at the same delayed coor- 
dinate itiib,x,z) = (fobs~ '^''f'' ,x,z) as the comoving density. 
Note that when off-axis calculations are performed (i.e., when 
the line of sight and the jet axis do not coincide), the photo- 
sphere location becomes a function of x and y and the above 
equation should be re-written with the explicit dependence on 
y as well as on x. 

Once the photosphere's coordinate Zph is obtained, we com- 
pute the light curves as (paper Q) 



Wobs) - ■ 



dxdy 



(l-/?cos6',)2 



(2) 



where T' = (3/?/fl)'/^ is the comoving temperature, a = 7.56 x 
10"'^ erg cm"^ K"'* is the radiation constant, and x^ax is a dis- 
tance perpendicular to the line of sight large enough so that 
the contribution to the emission atx^ax is negligible, y^ax sat- 
isfies the same constraint as JCmax- Finally, the peak frequency 
of the observed spectrum is computed as 

hu^,,,^ = 2.%6kBT' (3) 

were 5=1 /F(l - /3cosf?,,) is the Doppler factor. 

In principle, one would want the observer location Zobs to 
be approaching infinite, but practical considerations force it to 
be inside the computational area. This has the effect of bias- 
ing the photospheric distance, making it smaller than it should 
be. Assuming that the evolution at radii larger than the calcu- 
lated photosphere is self similar and that there is no further 
acceleration, one can easily evaluate the correction thanks to 
the fact that the opacity in a wind scales as the inverse square 
of the distance. Simple algebra yields: 



1 



■^ph.tme 
^ph,sim 1 ~ ■^ph.sim/'^obs 



(4) 



where Zph.tme is the true location of the photosphere and Zph sim 
is the photospheric radius evaluated from a simulation with 
observer location Zobs- Not surprisingly, the correction is 
small for Zph.sim < Zobs, but for Zph.sim = Zobs/2 the correc- 
tion is not negligible. Under the same assumption of non- 
accelerating self-similar expansion, a correction for the light 
curve and peak frequency can be obtained, taking into account 
that the photon-to-baryon ratio does not change with distance 
(e.g., Lazzati & Begelman 2010) and that the photon temper- 
ature scales as r~~l^ (e.g., Meszaros & Rees 2000). The true 
photospheric luminosity can be computed as: 



^true ~ ^sim (l Zph.sim/Zout) 



(5) 



Note that this equation is more precise than the one reported in paper 
I where the integration is en'oneously peii'ormed over an angular direction 
rather than a linear one. It is also generalized to include lack of any axial 
symmetry, as is the case when the jet axis and the line of sight do not coincide. 




t(s) 

Fig. 2. — Photospheric distance for the lines of sight analyzed in this paper. 
The photospheric distance is shown for the (x,y) = (0,0) location, i.e., along 
the line of sight. Since the distance of the virtual observer is Z^bs = 1.3 X 
cm, all the photospheric distances that ai'e measured to be close to .Zq^i; 
are likely affected by the finite size of the box (cfr. Eq.|4). 

and the photon frequency analogously scales as: 



t'true 



^1 Zph.sim/Zout) 



n2/3 



(6) 



Figure [T] shows the results of the light curve calculation at 
nine different observer angles 6*0 = 0, 1,2, 3, 4, 5, 6, 7, and 
8 degrees. The light curves at 6*,, = 2.5 and 3.5° were com- 
puted as well but are not shown in the figure. Light curves 
at 6'o > 8° were not computed because the photospheric ra- 
dius turned out to be too close to the outer boundary of the 
simulation box (see Figure |2] where all the photospheric radii 
are shown). Thick solid lines show the light curves as calcu- 
lated from the simulations (Eq. |2]i while thin lines show the 
light curves once the correction for the finite size of the box 
is applied (Eq.|5]). The comparison between the thick and thin 
lines shows that the correction is not large, mostly a factor 
2, but becomes progressively more important at large off-axis 
angles. The figure confirms the results of our paper I (see also 
Mizuta et al. 2010; Nagakura et al. 2010) that bright light 
curves can be produced by photospheric radiation. Compar- 
ing the light curves at different angles we see that there seem 
to be a sharp difference between the light curves at off-axis 
angles 6,, < 3 and those at 6'o > 3. The inner light curves 
have a bright early phase, characterized by some level of vari- 
ability, while the outer light curves have a dim early phase, 
becoming more luminous at times when the inner light curves 
have lost their brightness. The inner light curves are bright 
during the "shocked jet" phase, as discussed in Morsony et al. 
(2007), and become much more dim in the "unshocked jet" 
phase, during which they are flat and featureless (at f > 40 s). 
This difference is due to the fact that during the "shocked jet" 
phase the jet is hydrodynamically squeezed and its opening 
angle is smaller than the injection opening angle, so that very 
little emission is seen at > 9j ^ 4° (see Figure 6 of Morsony 
et al. 2010). We also note that the very sharp rise and decay 
times seen in the on-axis light curve are not present in the off- 
axis ones. We conclude that the 60 = 0° light curve is affected 
by 2D numerical artifacts and its temporal shape should not 
be taken as conclusive. We finally note that the sharp turn-off 
of the light curves at f ^ 95 s is due to the fact that we turned 
off the engine at f = 100 and is therefore an artifact of our 
choice for the engine duration. 

^ The ~ 5 s difference is the time is take the jet to break out of the progen- 
itor star. 
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Fig. 3. — Radiative efficiency of tlie pliotospliere's radiation as a function 
of tlie off-axis observer. Gray symbols sliow results extracted directly form 
the simulations, while black diamonds show the result when the coirection 
for the finite simulation box is applied. 

An important conclusion we reached in paper I was that 
of the large radiative efficiency of the photospheric emission. 
The radiative efficiency is computed as the ratio between the 
energy that is radiated at the photosphere and the total kinetic 
energy in the flow just below the photosphere. In this pa- 
per we confirm our previous finding, even though we point 
out that the efficiency computation is severely affected by 
the finite size of the computation box. Figure |3] shows the 
efficiency computed directly from the simulation (gray cir- 
cles) and the one derived by applying the finite box correction 
(Eq. |5]l. The efficiencies computed in the two ways are dif- 
ferent by a factor ^ 2. In any case, even the lower estimate 
puts the efficiency of the photospheric emission in the 30% 
range, larger than any efficiency that can be obtained in inter- 
nal shocks, without assuming an ad-hoc broad distribution of 
Lorentz factors (Lazzati et al. 1999; Panaitescu et al. 1999). 
The figure also shows the photospheric efficiency for a coni- 
cal jet with the same initial conditions that we adopted in our 
simulation (8%). Our simulation shows that the jet coupling 
with the progenitor star material has the effect of increasing 
the photospheric efficiency by a factor 3 to 8, with the high- 
est increase observed at relatively large off-axis angles (3 to 8 
degrees). The efficiency seems to have a somewhat bimodal 
behaviour, with the light curves for 9o < 4° having a ^ 30% 
efficiency and the light curves for 9o > 4° have a ^ 50% effi- 
ciency. The difference of light curves between small and large 
angles has been noted above. It seems that at small angles the 
light curve is dominated by the emission from the shocked-jet 
phase (Morsony et al. 2007), with an efficiency of ^ 30%. At 
larger off-axis angles, however, the emission is dominated by 
the phase during which the jet enters the corresponding line of 
sight. In this case the photospheric efficiency is larger, prob- 
ably as a cosnsequence of the dissipation due to the shearing 
motion of the jet with the progenitor star material. 

Another parameter we extracted from the simulations is the 
peak frequency of the light curves (Eq. O. Figure |4] shows 
the result of the temporal evolution of the peak frequency 
for the various lines of sight. First, we note that the ob- 
served peak frequencies are in fairly good agreement with the 
peak frequencies observed in GRBs, with an observed mode 
hi^obs = 250 keV (Kaneko et al. 2006). A moderate hard- 
to-soft evolution is observed for the 9,, = 0° curve. This is 
intriguing, given the so far unexplained hard-to-soft trend in 
many GRBs. However it could be an artifact of the 9o = 0° 
line of sight, since the trend is not observed at larger off-axis 
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Fig. 4. — Temporal evolution of the peak frequency of the light curves for 
the various off-axis angles. The 60 = 0° curve is shown only for the first 
33 seconds, since a sudden increase to a peak frequency of ~ 700 keV is 
seen afterwards. This sudden increase is due to the onset of the unshocked 
jet phase (Morsony et al. 2007), but the actual value of the peak frequency 
is a numerical artifact due to the coincidence of the line of sight with the 
symmetry axis of the simulation. 

angles. The comparison of the outer light curves (9„ > 3°) 
with their peak frequency evolution suggests the presence of 
tracking behaviour, i.e., increasing peak frequency when the 
light curve brightens and decreasing peak frequency when the 
light curve dims. This behavior is shown in Figure |5] for the 
9o = 4° light curve. 

A fundamental aspect of the observed ORB spectrum is 
the presence of non-thermal power-law tails at either end of 
the spectrum. The origin of the high-frequency tail in photo- 
spheric emission has been widely investigated and it has been 
shown that whenever some level of dissipation is present in 
the sub-photospheric region (10 < tj- < 1), inverse Compton 
processes produce a power-law tail comparable to the obser- 
vations, extending up to tens or hundreds of MeV (Pe'er et 
al. 2005, 2006; Giannios 2006; Giannios & Spruit 2007; Laz- 
zati & Begelman 2010; Beloborodov 2010). A more serious 
problem for photospheric emission models is the presence of 
non-thermal low-frequency tails. If the spectrum is gener- 
ated in a 1-zone photospheric model with radiation and mat- 
ter in thermal equilibrium, non-thermal low-frequency spec- 
tra are prohibited since they would violate the Raleigh- Jeans 
limit. Pe'er & Ryde (2010) explored the possibility of pro- 
ducing low-frequency non-thermal tails by considering the 
multi-color black body spectrum integrated over the emitting 
surface. They conclude that flat low-frequency spectra com- 
parable to those observed can be indeed produced by multi- 
color photospheres, but only during the decaying part of the 
light curve. Pe'er and Ryde (2010) considered emission from 
conical jets, for which analytic predictions can be made. We 
have computed the multi-color spectrum from our simulations 
to check whether a non-thermal spectrum at low-frequencies 
can be obtained during the raising part of the light curve, re- 
laxing the conical jet approximation. The resulting spectrum 
is shown in Figure|6] compared to a thermal spectrum with the 
same peak frequency (dashed line). As the figure shows, we 
do not find any significant deviation from a Planck law, be- 
cause the emission is dominated by the parcel of the outflow 
that moves towards the observer This conclusion is not in dis- 
agreement with the calculation of Pe'er & Ryde (2010), since 
they could find deviations from a thermal spectrum only when 
the material pointing directly at the observer is not producing 
any radiation. 
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Fig. 5. — Photospheric Hght curve for the 9a = 4° hne of sight (sohd Hne, 
left y axis) compared to the evolution of the peak frequency of the spectrum 
(dashed line, right y axis). 

4. A QUANTITATIVE COMPARISON 

In order to better evaluate the role of photospheres in the 
production of the prompt GRB radiation, we attempt a quan- 
titative comparison of our results with spectral and luminosity 
data. It is well known that observed GRBs for which a redshift 
measurement is available obey to the Amati correlation (Am- 
ati et al. 2002). According to this correlation, the burst frame 
peak frequency is correlated with the total isotropic equivalent 
luminosity elevated to an exponent close to 0.5. The most re- 
cent best fit to the correlation reads: 

(Amati, private communication), where the peak frequency is 
expressed in keV. The dispersion of the data around the best 
fit line is Gaussian with a sigma of 0.2. Figure [T] shows the 
best fit correlation with a dashed line and the 3 - u dispersion 
areas with dotted lines. The location of our light curves in the 
Amati plane is shown with gray circles (calculated directly 
from the simulation) and black diamonds (after the finite box 
correction). The simulation points are connected with a line 
showing the evolution from the on-axis light curve (upper 
right point) to the most off-axis light curve (lowest left point). 

The comparison of the simulated light curves with the Am- 
ati correlation is contradictory. On the bright side, as the ob- 
server moves away from the jet axis, the isotropic energy and 
peak frequency evolve in such a way as to reproduce the slope 
of the correlation. This is independent of whether we consider 
the corrected or uncorrected light curves as the most reliable 
result. On the other hand, the normalization is incorrect, since 
the simulated light curves are too bright for the corresponding 
peak frequency (see also Nagakura et al. 2010). 

5. DISCUSSION 

This paper presents the natural extension of our 2009 
numerical study of the photospheric emission from long- 
duration GRB jets (paper I). We here present a simulation run 
on a domain ten times bigger than in paper I. Such a simu- 
lation allows us to explore the location, spectrum, and light 
curves of photospheres along lines of sight different from the 
jet axis. The results of this paper substantially confirm the 
results of paper I, even though we find that some details of 
paper I had been affected by the relatively small computa- 
tional domain and by the coincidence of the line of sight with 
the symmetry axis of the computation. For example, the sharp 



Fig. 6. — Observed spectrum for the 6„ = 1° observer at f = 30 s, integrated 
over the emitting surface (upper panel, thick solid line). The spectrum shows 
only small deviations with respect to a Planck function with the same maxi- 
mum (thin dashed line). The bottom panel shows the local power-law slope 
of the two spectra (L(v) oc i/°), emphasizing that the difference in the slopes 
is very small. 

features seen in the on-axis light curve disappear as soon as 
the observer moves away from the jet axis. However, we find 
that the high efficiency of the photospheric radiation is con- 
firmed and that the efficiency increases for off-axis lines of 
sight. 

Even with the big simulation box presented here (Zom = 
2.5 X lO'^ cm), the photosphere during some times interval 
is sometimes affected by the box size. We therefore discuss 
a correction that allows us to derive the correct location of 
the photosphere under the approximation of self-similar non- 
accelerating expansion beyond the box outer edge. Off-axis 
photospheric light curves were discussed elsewhere (Mizuta 
et al. 2010; Nagakura et al. 2010). Mizuta et al. perform a 
simulation with an identical setup to our paper I. They derive 
off-axis light curves and spectra, but their simulation box has 
an outer edge Zout = 2.5 x 10'^ cm, making their results unre- 
liable. In fact, our Figure|2]shows that for all lines of sight the 
photosphere is well outside their simulation box. Nagakura 
et al. (2010) perform instead a somewhat low -resolution sim- 
ulation extending to very large radii. While their simulation 
is superior to the one we present here in terms of size, their 
resolution and the adoption of a fixed-grid algorithm make 
their results inaccurate in terms of the internal properties of 
the outflow. In particular, a fixed-grid algorithm is not able to 
resolve the tangential shocking that we observe in the jet out 
to hundreds of stellar radii, losing the dissipation due to such 
shocks and therefore underestimating the temperature at the 
photosphere. 

The size of this simulation box allows us to perform an in- 
teresting comparison between the physically motivated pho- 
tospheric light curves and the light-power curves that we have 
used in other publications as proxies to the real light curve 
(Lazzati et al. 2010ab; Morsony et al. 2010). The comparison 
is shown in Figure[8]for the 6*0 = 1° line of sight. Light-power 
curves were extracted at radii of 10' ' and 10'^ cm, assuming 
an efficiency of 30%. We find that in all cases the three bolo- 
metric light curves show a strong correlation, demonstrating 
that the use of light-power curves is justified. However, the 
light-power curve extracted at the largest radius does show 
a bright early spike that is not present in the other curves. 
This spike is due to the presence of slow, dense material that 
piles up ahead of the jet and it may well be an artifact of the 
presence of a reflecting boundary condition along the jet axis. 
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Fig. 7. — Location of our simulated photospheric light curves and spectra in 
the Amati diagram. The dashed line shows the location of the best-fit Amati 
relation, while the dotted lines show the 3 — cr dispersion ai'ound it. Black 
and gray symbols show the results with and without the finite-box correction, 
respectively. 

The comparison of photospheric hght curves with light-power 
curves at different off-axis angles produces results in qualita- 
tive agreement with what shown for the 9o= 1° case. 

Two of our new and more interesting results (not discussed 
in paper I) concern the spectrum of the photospheric emission. 
First, we find that the integration of the multi-color black- 
body spectrum over the emitting surface does not alter sub- 
stantially the Planck law shape of the photospheric emission 
(Figure |6]l. Second, we find that the peak frequencies pro- 
duced by the photosphere are too small for the isotropic equiv- 
alent energy of their respective light curves (Figure|7]i. These 
two results are important because they show that the photo- 
spheric emission has trouble in reproducing the non-thermal 
low-frequency spectra observed in GRBs and their peak fre- 
quencies. Both these problems may be solved if the assump- 
tion of thermal equilibrium between the baryons and the radi- 
ation field that we adopted to derive the spectrum are inade- 
quate. As a matter of fact, the bolometric light curve calcula- 
tion relies on the simple assumption that the internal pressure 
of the outflow is mainly due to radiation and is therefore quite 
robust. Such assumption is correct no matter the spectrum 
of the radiation field. If, however, the radiation has a color 
temperature that is higher than its effective temperature, our 
spectral calculations would be incorrect. A first consequence 
would be that the peak frequency that we derive is underesti- 
mated, and therefore the Amati correlation discrepancy could 
disappear As a matter of fact, in scattering dominated accre- 
tion disk atmospheres, a discrepacy of a factor ~ 3 is found 



in detailed radiation transfer calculations (e.g., Davis et al. 
2005). In addition, having a color temperature higher than the 
effective temperature, allows for the presence of non-thermal 
low-frequency tails without violating the Raleigh- Jeans limit. 

Even though finding the exact mechanism that would cre- 
ate such an off-balance configuration between the baryon and 
the radiation field is beyond the scope of this paper, we argue 
that research on the photospheric emission of GRB jets should 
proceed in two directions. On the one hand, it is mandatory to 
explore the shape of the spectrum of the photospheric emis- 
sion from first principles, relaxing the assumption of equilib- 
rium. On the other hand, the role of different jet and progeni- 
tor properties should be explored, since all the available simu- 
lations have so far concentrated on the progenitor model 16TI 
(Woosley & Heger 2006) and on jet configurations analogous 
to the one we exnlored in naner I. 
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Fig. 8. — Comparison between the 6o = 1° photospheric light curve (thick 
solid line) and light-power curves extracted at the same viewing angle at radii 
r = lO" cm (thick dashe fine) and r = lO" cm (thin solid line). 
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